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I ^ I Abstract 

^ We report the final result of the CUORICINO experiment. Operated between 2003 and 2008, with a total exposure of 19.75 kg-y 
j> of '^"Te, CUORICINO was able to set a lower bound on the '^"Te Ov/3/3 half-life of 2.8 x lO^^ years at 90% C.L. The limit here 
\^ reported includes the effects of systematic uncertainties that are examined in detail in the paper The corresponding upper bound on 
the neutrino Majorana mass is in the range 300-710 meV, depending on the adopted nuclear matrix element evaluation. 

(N 

cn 

CN 1. Introduction ery potential [1, 2]. 
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The Ovf3/3 transition could proceed through several different 
Double beta decay (J3/3) is a rare transition between two iso- mechanisms, of which the simplest and most commonly cited 



• • bars, involving a change of the nuclear charge Z by two units. is the exchange of a light Majorana neutiino. In this case the 

. ^ There exist several naturally-occurting even-even nuclei for observation of Ov/3/3 would not only provide evidence for lep- 

^ which this is the only allowed decay mode. While the transition ton number violation and the Majorana character of the neu- 

involving 2 electrons and 2 (anti)neutrinos (2v/3/3) is allowed in trino, but would also result in a measurement of the effective 

^ any theoretical model (and has been obsei-ved for various iso- Majorana mass - I 2 miU^.l (where m,- are the three mass 

topes), this is not the case for the neutrinoless channel {Ov/3/3). eigenstates of neutrinos and Uei are the PMNS matrix elements) 

Despite being energetically possible, Ov/?/? violates lepton num- through the relation: 
ber conservation by 2 units and is possible only if the neutrino 

is a massive Majorana particle [1]. Ov/3/3 searches have been 1 _ ^2 _ ^2 q'-^v^j^Ov^i 
pursued for more than half a century and today they experience Tj '''' 
a renewed interest, thanks to the discovery of neutrino oscilla- 
tions [2, 3, 4]. Neutrino oscillations imply that neutrinos have Here t°^2 the Ovfi/3 half-hfe, G"" is the two-body phase-space 
a finite mass, but are not enough to determine the nature of that factor and M°'' is the Ov/3/i nuclear matrix element (NME). The 
mass (Dirac or Majorana) or why it is so extraordinarily small. product = G*'''|M"''|^ is called the "nuclear factor of merit"; 
Several theoretical speculations point toward a mass generation the name refers to the fact that, according to equation ( 1), di- 
mechanism that implies a Majorana character of neutrinos, and rectly influences the experimental sensitivity to m^e- The main 
that indicates the OvySyS process as the unique tool with a discov- uncertainty in deriving mj,? (or an upper limit on it) from the 
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Figure 1: Anticoincidence total energy spectrum of all CUORICINO detectors (black). The most prominent peaks are labeled and come from known radioactive 
sources such as: e+e" annihilation (1), -'''Bi (2), **''K (3), '"''TI (4), "^"Co (5) and Ac (6). The total energy spectrum of all CUORICINO detectors during 
calibration measurements is also shown (color). For convenience, it is normalized to have the same intensity of the 2615 keV line of 20871 as measured in the 
non-calibration spectrum. 



experimental result on r^J^2 comes from the NME, which is a 
theoretical calculation still affected by a large spread among the 
adopted nuclear models and their implementations [5, 6, 7, 8]. 
To mitigate this uncertainty, several candidate Qvfifi isotopes 
could be studied and the m^^ results compared. Although it 
may not be feasible to study all the /3p candidates with reason- 
able sensitivity, there exists a "golden list" of isotopes that - as 
a compromise of cost, availability, technological approach and 
other factors - have been studied so far Among them, those that 
have yielded the most stringent limits on niee (within the NME 
spread) are ^^Ge [9, 10], i°"Mo [11], '^"Te [12] and '^^Xe [13]. 
In all but one [14] case, only upper bounds on the Majorana 
mass have been reported. 

In this paper, we discuss the final Ov/3/3 result of CUORI- 
CINO, which yields the most stringent bound on based on 
'^"Te studies, and one of the best in general. CUORICINO data 
acquisition started in April 2003 and ended in June 2008. The 
data are separated into two runs (RUN I and RUN II), due to a 
major maintenance interruption. The data collection is summa- 
rized in table 1. A partial data-set of 11.83 kg-y of '^"^Te expo- 
sure was used for the analysis reported in [12]. 

The paper is organized as follows: after a short description 
of the experimental set-up in section 2, we present details of 
RUN II data analysis, discussing processing in section 3, data 
reduction in section 4 and efficiency evaluation in section 5. 
In section 6, we describe two Bayesian approaches used for 
the Qv/3/3 half-life limit evaluation (one of them is that adopted 
in [12]), testing the two procedures on a toy Monte Carlo and 
discussing their compatibility on real data. In section 7 we dis- 
cuss the influence of systematic errors on the final result. 

We conclude the paper with the CUORICINO final result for 
the OvjSjS half-life limit evaluated on the entire data set. This 
is done treating RUN I and RUN II as two independent ex- 
periments whose likelihoods are combined. This choice was 
motivated by the difference in detector configuration between 
the two runs (increased number of active detectors, improved 



performance) and a presumable difference in background com- 
position (due to detector exposure to air). 



2. CUORICINO 

For many years, the most sensitive Ov/3/3 results for '^"Te have 
been obtained using bolometric detectors. A bolometer is a 
type of calorimeter operated at ultra-low temperature, in which 
the energy of incident radiation is converted to heat, raising 
the temperature of the detector's body. The energy released 
in the detector is then determined by measuring the tempera- 
ture increase, in our case by using a semiconductor thermistor 
[16] whose resistance varies exponentially with temperature. In 
this kind of detector, the candidate isotope is contained within 
the active mass of the detector itself. CUORICINO bolometers 
contain the isotope '^"Te (isotopic abundance 33.8%) which is 
a /3/3 candidate with a rather favorable factor of merit' . 

About 85% of the time (see section 4), the two electrons 
emitted by '^°Te Qv/3/3 would be fully contained within one 
crystal. The signature of the decay would therefore con- 
sist of a monochromatic peak in the energy spectra of the 
bolometers at an energy equal to the Q-value of the decay: 
2527.518±0.013 keV [17]. The difficulty of the experiment 
lies in the control and reduction of all the background events 
that could mimic such a signal. These can be non-particle sig- 
nals, due to electronic or thermal noise, or particle signals, due 
to radioactivity and cosmic rays. The former are rejected on 
the basis of a pulse shape discrimination technique (see sec- 
tion 4). The latter are controlled during the experiment's de- 
sign and construction by proper material selection, handling, 



'a second fSfi candidate with high natural abundance is the isotope '^^Te 
(isotopic abundance 31.7%). This isotope is not as interesting as '^"Te because 
of its lower transition energy, which reduces the nuclear factor of merit and 
shifts the signal to a higher background region (therefore also lowering the 
achievable experimental sensitivity). 
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Table 1: CUORICINO crystal information and statistics. Crystal mass is the average measured mass for CUORICINO detectors. 



Crystal Type 


Crystal Mass 
[g] 


i3"Te Mass 
[g] 


Exposure Run 11 
[kg(i30Te)-y] 


Exposure Run I 
[kg(i30Te).y] 


big 


790 


217 


15.80 


0.94 


small 


330 


91 


2.02 


0.094 


'^"Te-enriched 


330 


199 


0.75 


0.145 



and shielding [12, 18, 19], and — at the stage of data-analysis 
— by coincidence cuts (see section 4). 

CUORICINO was the latest step in a long series of exper- 
iments of increasing mass, performance and sensitivity [20, 
21, 22]. The next experiment to use this technique will be 
CUORE [23], which is presently under construction. 

CUORICINO [12] was a tower array of 62 TeOi crystals 
used as bolometric detectors. The array was operated un- 
derground, in a dilution refrigerator located at the Laboratori 
Nazionali del Gran Sasso (INFN - Italy), which provides an av- 
erage coverage of 1400 m of rock (3650 m.w.e.). CUORICINO 
crystals can be divided into four main groups according to their 
mass and isotopic abundance. These are: 

• the "big crystals" — 44 bolometers, 5x5x5 cm^ in size 
and 790 g in mass; 

• the "small crystals" — 14 bolometers, 3x3x6 cm-' in size 
and 330 g in mass; 

• the "'^"Te-enriched crystals" — 2bolometers, 3x3x6cm^ 
in size and 330 g in mass, grown with '^"Te enriched ma- 
terial; 

• the "'^**Te-enriched crystals" — 2 bolometers, with the 
same size and mass of the '^"Te-enriched ones but grown 
with '^^Te enriched material. 

The '^"Te-enriched crystals have a '"'"Te content corresponding 
to an isotopic abundance of 75% [22], while the '^"Te content of 
128-enriched crystals is so low that they will not be considered 
for the OvjSjS analysis presented in this paper 

A detailed description of the array, cryogenic set-up, shields, 
front-end electronics, and DAQ can be found in [12] and ref- 
erences therein. Here we will describe the main steps of data 
acquisition and data handling which are relevant for the discus- 
sion. 

The output voltage of each detector was monitored by a con- 
stant fraction trigger When the output voltage exceeded the 
trigger threshold, the acquisition system recorded 512 sam- 
ples (a 4 s window sampled at 125 Hz), which constitute one 
"event". The acquired time window fully contained the pulse 
development, providing an accurate description of its wave- 
form. A pre-trigger interval just prior to the onset of the pulse 
was used to measure the DC level of the detector, which cor- 
responds to the instantaneous detector temperature. The pulse 
amplitudes were evaluated offline for each recorded event, to- 
gether with a few other characterizing parameters of the pulse. 
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Figure 2: A bolometric particle event (left) and a spurious non-physical signal 
(right), superimposed with their optimum filter output (in color) in the time 
domain. 



Each single CUORICINO measurement lasted about 22 
hours on average, with the time between measurements (about 
2 hours) dedicated to cryogenic system maintenance. A routine 
calibration with an external ^^^Th source was performed 
approximately once per month, lasting for about 3 days. The 
accumulated data between two calibrations is referred to as 
a "data-set". The spectrum obtained by summing all the 
CUORICINO collected data (i.e. summing over detectors and 
data-sets) is shown in figure 1. The background recorded by 
the detectors is clearly dominated, in this region, by gamma 
emissions due to radioactive contaminations of the detector 
and of the surrounding apparatus. The most intense gamma 
lines are listed in reference [12]. In figure 1, the spectrum 
corresponding to the sum of all calibration data is also shown. 
For convenience the calibration spectrum is normalized to have 
the same intensity of the 2615 keV line of -'^'^Tl as measured in 
the background spectrum. 



3. Data processing 

The analysis of CUORICINO data starts with the collection 
of all the triggered events. For clarity, we will model the single 
waveform V{t) induced by a particle interaction in the crystal 
as: 



V(t) 



Vo s(t) + n(t) 
G{T)-A(E) 



(2) 
(3) 



In the first equation, Vq is the maximum value of the raw sig- 
nal acquired at time fo, s{t) describes the shape of the particle 
signal and n(f) is an additive noise source. The second equation 
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describes the dependence of the signal ampHtude on the detec- 
tor working temperature. Here we assume that the dependence 
of the gain on temperature, G(T), and that of the ampHtude on 
energy, A(E), can be factorized. This is not true in general, how- 
ever it is a good approximation when dealing — as in our case 
— with small temperature drifts. Although it describes a naive 
model, this formula highlights the key points of the analysis. In 
order to estimate E, we need: 

1. a technique able to measure Vq, reducing the effect of «(f) 
as much as possible, in order to improve our resolution 
(amplitude evaluation); 

2. an algorithm to control for the variation of G{T) produced 
by detector temperature drifts (gain instability correction); 

3. a technique to measure the form of A(E) (energy calibra- 
tion). 

Amplitude evaluation. This is done by maximizing the 
signal-to-noise ratio by means of optimum filtering [24]: each 
waveform is convolved with a transfer function li(t) whose 
Fourier transform is defined as: 



H((^) = e" 



S*iaj) 



(4) 



N(aj) 

where S (w) is the Fourier transform of the average detector re- 
sponse function s{t), N((d) is the spectral power density of the 
noise characterizing the detector, and t,„a^- is the time at which 
the pulse reaches its maximum. The functions s{t) and N(a)) are 
computed by an averaging procedure of the bolometric pulses 
and of the Fourier transformed baselines^. 

Figure 2 shows an example of an event due to a particle 
interaction and and example of a non-particle event, most likely 
due to an abrupt temperature increase produced by an electric 
disturbance or by vibrations. Each waveform is superimposed 
with its optimum-filtered counterpart. 

Once the optimum filter is applied, the amplitude of the signal 
is inferred from the maximum of the filtered waveform in time 
domain. 

Gain instability correction. This correction is achieved by 
measuring the voltage amplitude, Vref, of a monochromatic 
reference pulse. This pulse is produced by depositing a fixed 
amount of energy into the crystal by the Joule dissipation from 
a heavily doped silicon resistor glued to the crystal. Because 
the energy deposited is fixed, any variation of Vref would be 
due to a variation in G{T), which can be measured and used to 
correct the amplitudes of all the triggered events. For a more 
detailed discussion of this method we refer to [25, 26]. 

Energy calibration. The voltage-energy relationship is re- 
constructed by means of routine source calibrations: two wires 
of thoriated tungsten are periodically inserted between the cryo- 
stat and its external lead shield. The voltage amplitude of the 
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Figure 3: Residuals (nominal energy - calibrated energy) vs. nominal energy 
evaluated on the main gamma lines identified in the CUORICINO background 
spectrum. Circles (in color) refer to a calibrated energy obtained with the third 
order polynomial, triangles (in black) with the log-polynomial. 



pulses corresponding to the main gamma lines of ^^^Th are 
used for the determination of the parameters describing the 
A(E) relationship. This function is characterized by different 
non-linearity sources [12], the dominant one being the depen- 
dence of the thermistor resistance on the temperature [27]. In 
this work, A{E) is parametrized with a third-order polynomial, 
which can be considered as the truncated Taylor's expansion of 
the real unknown calibration function. In the previous CUORI- 
CINO analysis [12], a different calibration function was used. 
This was a second order polynomial in log(y) and log(£'), based 
on a thermal model describing our bolometric detectors. While 
this function performs better at extrapolation (i.e. above the 
highest calibration line at 2615 keV), the third order polyno- 
mial performs better in the interpolation region (i.e. between 
threshold and 2.6 MeV). The difference between the two pa- 
rameterizations was studied using the total background spec- 
trum recorded by CUORICINO (figure 1); this spectrum con- 
tains several gamma peaks whose origins, and therefore nomi- 
nal energies, are clearly identified. The difference between the 
nominal energy of each peak and its measured position (the 
residual) is plotted against the nominal energy in figure 3, show- 
ing the slightly better performance of the 3rd order polynomial. 
These residuals also provide important information concerning 
the precision of our calibration: their spread can be used as an 
estimator of the uncertainty in the energy position of a peak, 
including that produced by the Qv/3/3 signal. 

Source calibration measurements are repeated for each data- 
set and are also used to check the detector performances over 
time. Figure 4 shows the distribution of all the resolutions mea- 
sured in calibration for the three crystal groups (big, small, and 
enriched crystals). 



-At random times during the course of data acquisition, sets of 5 12 samples 
("baselines") were collected in anticoincidence with the trigger. These were 
used for the evaluation of the average noise power spectrum, as described in the 
text. 



4. Data reduction 

The final CUORICINO spectrum is composed of events 
which survived two different types of data selection: global 
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Figure 4: Distribution of the energy resolutions (FWHM) measured in calibra- 
tion for the three groups of crystals during the 33 data-sets belonging to RUN 
II. From top to bottom: big crystals, small natural crystals and '^"Te enriched 
crystals. 

and event-based cuts. 

Global cuts: these are applied following quality criteria 
decided a priori (e.g. an excessive noise level or an incom- 
patibility between the two calibration measurements at the 
beginning and the end of a data-set). They identify bad time 
intervals to be discarded. This kind of cut introduces a dead 
time that is accounted for by properly reducing the live time 
of the detector of interest. The cuts are generally based on 
off-line checks that monitor the detector performances and flag 
excessive deviations from global control quantities (average 
resolution, average rate, etc.). The total dead time introduced 
by these global cuts is ~5%. A further dead time is introduced 
by the rejection of a short time window centered around each 
reference pulse (the frequency with which the reference pulses 
are generated is about 3 mHz). This cut ensures the rejection 
of possible pile-up of a particle signal with the reference pulse 
(the impact of this cut is reported as an efficiency in table 2). 



Event-based cuts: these are the pulse-shape and the anti- 
coincidence cuts. The former is used to reject non-physical 
and pile-up events (the presence of a pile-up prevents the op- 
timum filter algorithm from providing a correct evaluation of 
the pulse amplitude). The latter allows for the reduction of the 
background counting rate in the region of interest (ROI). The 
Ovy6/3 signature we look for consists of a "single-hit" event (only 
one detector at a time involved), while many of the background 
counts in the ROI are due to "multiple-hit" events. These in- 
clude events due to alpha decays on the crystal surfaces that 
deposit energy in two neighboring crystals and events due to 
gammas that Compton scatter in one crystal before interacting 
in another one. 

The pulse shape parameters used in this analysis are the rise 
time and decay time of the raw waveform and a parameter that 
measures the consistency of one of the basic statements of opti- 
mum filter theory. This "Optimum Filter Test" parameter (OFT) 
is the difference (expressed in percentage of the total amplitude) 
between the evaluation of the pulse height in the time domain 
(as the maximum value of the filtered pulse) and that in the fre- 
quency domain (as the integral of the filtered-pulse power spec- 
trum). Indeed, if the shape of an event is identical to the average 
detector response, the two methods yield the same result. How- 
ever, if the shape of the signal is different from the expected 
one (such as in the case of a non-physical or a pile-up event), 
they differ. Figure 5 shows the scatter plot of OFT as a function 
of energy for a CUORICINO detector (here only one data-set 
is reported). The main trend reflects the change of the signal 
shape with energy, and has a minimum in the region where the 
average response was measured (1-2 MeV). This variation in 
the pulse shape is caused by non-linearities introduced by the 
thermistor Outliers on this plot correspond to misshapen events 
which will be discarded; the colored vertical bars identify con- 
fidence regions evaluated automatically on this distribution in 
order to obtain cuts which are independent of the signal ampli- 
tude. 

The anticoincidence cuts require that only one detector fires 
within a time window of 100 ms. 

5. Signal efficiency 

The signal efficiency is the probabihty that a Qv/3/3 event 
is detected, its energy is reconstructed accurately, and that it 
passes the data selection cuts. This parameter must be ac- 
curately determined, since it is used to obtain the number of 
^^°Te Qv/3/3 events. The overall signal efficiency is: (82.8 + 
1.1)% for the CUORICINO big crystals and (79.7 + 1.4)% for 
the small and the '^'^Te-enriched ones. These efficiencies were 
computed as discussed below. 

There are two main sources of inefficiencies, one "physical" 
that can be computed by simulations, and the other "instru- 
mental" that must be measured from the data. The mechanism 
of "physical" efficiency loss is the escape of a fraction of 
the OvySyS energy from the source crystal. Mechanisms for 
the "instrumental" efficiency loss are: the pulse-shape cut, 
the anti-coincidence cut and an incorrect assignment of the 
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Figure 5 : Typical scatter plot of the OFT (deviation of filtered raw signal from 
the average detector response) as a function of the signal's evaluated energy. 
The main trend identifies "good" events. Pile-up and non-physical pulses are 
outside the confidence regions (identified by the colored vertical bars). At low 
energy, there is a high density distribution of events for OFT values between 
10"" and 10"' : these events are due to electric disturbances. 



Table 2: Contributions to the CUORICINO Ov/3/3 signal efficiency. 



Source 



Signal efficiency (%) 



Energy escape 

Pulse-shape cuts 
Anti -coincidence cut 
Noise 



87.4+1.1 (big crystals) 
84.2+1.4 (small crystals) 

98.5+0.3 
99.3+0.1 
99.1±0.1 



Pile-up with reference pulses 97.7 



Total 



82.8±1.1 (big crystals) 
79.7±1.4 (small crystals) 



photopeak at 2615 keV due to ^osjj .^^^ ^^g^j ^ proxy for 
the Qvpp peak. The 2615 keV peak was chosen because of its 
proximity to the Qv^p energy and its relatively high intensity. In 
principle, an ideal pulse-shape cut should leave the main peak 
untouched and should only reduce the flat background. The 
area of the peak can then be computed in terms of the the total 
number of signal events (Nsig), the signal efficiency (eps), the 
total number of background events (Ntkg) and the background 
efficiency (ehkg)- A simultaneous fit was done on both the spec- 
tra of accepted and rejected events. The area of the peak in the 
accepted events spectrum is given by eps x Nsig, while for the 
rejected events it is (1 - eps) x A^,,,,. Similarly, the background 
yield for the accepted events is x Nbkg, and the background 
yield for the rejected events is (1 - £hkg) x Nbkg. By including 
eps directly in the parametrization of the fit, correlations among 
the fit parameters are automatically taken into account when the 
error on £ps is calculated. The result is eps - (98.5 + 0.3)%, 
and €hkg - (64 ± 2)%. The pulse shape cut is clearly very pow- 
erful, rejecting approximately 36% of the events in the contin- 
uum background while retaining 98.5% of the signal events in 
the peak. The events discarded in this region are mainly pile-up 
with real or spurious signals. 

To estimate the efficiency of the anti-coincidence cut, we 
used the same procedure but considered the only available high 
intensity peak that is produced by a nuclear decay with no de- 
tectable coincidence radiation. This is the 1460 keV gamma 
line emitted in "^''K electron capture (the only coincident radia- 
tion, a 3 keV X-ray, is far below our threshold). Since events 
in this photopeak are single hits, their reduction after an anti- 
coincidence cut can be ascribed only to random coincidences. 

The last source of inefficiency is the loss of Qv/3/3 events due 
to excess noise which can distort the pulse shape and introduce 
an error in the reconstructed energy. If such an error is greater 
than the resolution, the event can be considered as lost in the 
continuum background. In order to estimate this efficiency we 
compare the number of reference pulses generated during the 
measurements (the signals used for the gain instability correc- 
tion, see section 3) with the number actually measured in the 
correct energy range. 



energy of the signal (mainly due to noise and pile-up). Their 
contributions to the total signal efficiency are summarized in 
table 2. 

Physical efficiency: the OvySyS signature is a sharp peak cen- 
tered at the transition energy (Q-value) of the decay. The peak 
is produced by Qv/3/3 decays fully contained within the source 
crystal. The containment probability was evaluated using a 
Geant4-based Monte Carlo simulation that takes into account 
all the possible energy escape mechanisms (i.e. electrons, X- 
rays or bremsstrahlung photons escaping from the source crys- 
tals). Since the escape probability depends on the crystal geom- 
etry, the efficiency is slightly different for the big and the small 
crystals (see table 2). 

Instrumental efficiency: this is the product of the pulse- 
shape cut, anti-coincidence and excess noise efficiencies. To 
evaluate the efficiency of the pulse shape cut, the background 



6. analysis 

The definition of the energy window used to fit the Qvf)/3 spec- 
trum, the hypothesis assumed for the background shape and the 
number of free parameters used to describe the background it- 
self are extremely important for the choice of the analysis pro- 
cedure and for the determination of its systematic s. The choice 
of the energy window is somewhat arbitrary, but it influences 
the background representation. If the energy window is too 
wide (compared to the signal FWHM) a very precise knowl- 
edge of the background shape is necessary. Obviously there 
is also a minimum width necessary to be able to evaluate the 
background level beyond the Qv/3/i peak. In our case, there is 
a background line near the Ov/3/3 energy, at ~2505 keV, due to 
^"^Co (sum of the two photons emitted in cascade by ^''Co de- 
cay), which should also be included in the window. Given these 
considerations, our final choice for the fit window is 2474—2580 
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Figure 6: A closer view of the CUORICINO anticoincidence spectrum (pre- 
sented in figure 1 ) near tlie Oy/3/J ROI. Note tliat tlie efficiencies listed in table 2 
are not yet included. 



keV. This is the widest window centered on the Q-value that 
allows the following two background peaks to be excluded from 
the fit: the 2448 keV line of ^i^Bi and the 2587 keV Te X-ray 
escape peak of the ^^^Tl line. The latter peak is clearly visible 
in the CUORICINO calibration spectrum shown in figure 6, al- 
though — due to the lower statistics — it is not visible in the 
background spectrum shown in the same figure. 

Within this picture, we choose the simplest possible model 
for the shape fijiE) of the spectrum (normalized to mass, live 
time, efliciency and isotopic abundance) of each single detector 
(index /) in each data-set (index j) as: 



fij(E) = B,- + r^J" g^j {e - -t- gij {e - (5) 

Here gij(E) is the function describing the shape of monochro- 
matic energy lines in the detector, during the data-set, i.e. 
the response function that is represented by a gaussian with a 
width determined from calibration data^. E'^^ is the '-"'Te /3/3 Q- 
value, fixed at its measured value (2527.5 keV). E"-" is the sum 
energy of the two ''°Co gamma lines (2505.7 keV). B, is the flat 
background component for the detector (here we assume a 
time independent background). Finally T. j^" and r"''^'' are, re- 
spectively, the ^'^Co activity for the detector during the 
data-set (^"Co has a half-life of 5.27 years), and the absolute 
activity for Qv/3/3 , both expressed in counts/kg/y. 

Free parameters are B,, F^. F*''' and E"-", the parameters of 
the response function being fixed at the values measured during 
calibrations. Note that the dependence of F. on the index j is 
determined by the ^"Co half-life; therefore, the total number of 
free parameters is determined only by the number of detectors 
(i.e. by the index /)■ 
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Figure 7: Results of a toy Monte Carlo simulation with no Ov/?yS signal (1000 
simulated CUORICINO-RUN II experiments). Scatter plot of the 90% C.L. 
limits (top panel) and of the best fits (bottom panel) obtained with the two 
different approaches. The colored line has slope=I and shows the strong corxe- 
lation between the two techniques. 



6.1. Statistical approaches 

The CUORICINO spectrum shows no evidence of a 
OvySjS signal, thus we will provide a limit for the half life of 
^^°Te by means of a Bayesian approach. Unlike other low statis- 
tics methods such as that of Feldman and Cousins [28], this 
technique does not require an exact evaluation of the expected 
number of background events (which is unknown). In our case, 
all the uncertainties are marginalized in the process of the limit 
computation, and our prior knowledge for the rate F"'' will be 
represented by a flat distribution, excluding the non-physical 
region. 

Once the statistical method is chosen, we need to decide 
how to model the experiment: every CUORICINO detector 
can in fact be imagined as an independent search for OvySyS, 
characterized by its own background and resolution. There are 
three natural approaches which can be chosen to search for a 
0vy6/? signal: 



'To evaluate the energy resolution in the Oy/?/J region we use the 2615 keV 
peak since this is the nearest peak to the Q-value clearly visible in our calibra- 
tion spectra. 



I — treat the different detectors separately; 

II — treat the diff'erent detectors separately, assuming an 
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Figure 8: Results of a toy Monte Carlo simulation with no Ov/?/? signal (approach II in black and approach III in color). Left panel: pull distributions of the obtained 
best fits divided by their statistical eiTor. Right panel: distribution of the 90% confidence level limit on the decay rate. 



identical background for all detectors within each group 
(big, small, '^"Te-enriched); 

III — sum the spectra of all detectors belonging to the same 
group. 

In approach I, each detector and each data-set is fit with its own 
function fij(E). In principle this is the best approach since it 
uses all the information available; however, the number of free 
parameters is huge (about a hundred). 

Approach II lowers the number of free parameters by forcing 
the background and the ^'^Co rates to be identical on detectors 
of the same group. In this approach each detector and data-set 
is still described individually by its own function fijiE) but the 
total number of free parameters is reduced since B,- can assume 
only 3 values (for big, small and '^°Te-enriched crystals) and 
the same is true for r'j^"- This could be considered a strong 
assumption, but it is motivated by the fact that the low statistics 
prevent us from being sensitive to background variations among 
crystals of the same group in the Qv/3/3 region. This method also 
offers the advantage of being less sensitive to fluctuations in 
the counting rate of a single detector over time, and takes into 
account the decay rate of ^'^Co. 

Approach III removes the background assumption of the pre- 
vious model, at the price of a certain degree of information 
loss. The counting rate is simply averaged over all data-sets 
and detectors for the three mentioned groups. A variation of 
the background and of the ^"Co rate over time is then irrele- 
vant, provided that the response function does not change with 
time. The average is done simply by summing over all the data 
collected with detectors belonging to the same group, thereby 
obtaining three spectra that can be represented by the function: 

ME) = + r^C" G, {e - E^") + r"" G, {e - E°^) (6) 

Here, the index k has three allowed values for big, small and 
'^°Te-enriched crystals while the response function GkiE) is de- 



Table 3: Comparison between the fit results for the two studied approaches on 
RUN II data. 





Method II 


Method III 


Best Fit [y-'] 


(0.2+ 1.5)xlO-25 


(0.3 + 1.5)xlO-25 


Half-life limit [y] 


2.5 X 10^4 


2.4 X 10^4 



fined as: 

where the sum over / extends on all the detectors belonging to 
the k''' group and j runs over all the data-sets. Aij and cr,j are 
the corresponding background exposure and energy resolution 
measured during calibration. Note that — as is true also for 
the other two approaches — the response function is built using 
measured quantities, i.e. it does not contain any free parame- 
ters. 

We discarded the first approach due to the excessive number 
of free parameters, and we performed two parallel (and inde- 
pendent) analyses on real and Monte Carlo-simulated data for 
the limit computation, following an unbinned likelihood tech- 
nique [29] for the second approach and the standard CUORI- 
CINO Likelihood-Chi-Square technique [12, 30] for the third 
one. The goal was to choose the most reliable procedure, check- 
ing for possible biases and comparing performances. 

6.2. Comparison of methods on real data and Monte Carlo 
simulations 

Applying the two methods on the CUORICINO RUN II data, 
we obtained results which were indeed very similar. Table 3 
shows the best fits and limits for the two approaches: no signif- 
icant difference is observed. 
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Figure 9: Left panel: best fit, 68% and 90% confidence intervals for the total statistics (RUN I+RUN II) superimposed on the CUORICINO sum spectrum of the 
three groups of crystals (each scaled by efficiency and exposure) in the Ov/?/? region. (The purpose of the plot is to give a pictorial view of the result; in fact the fit 
was performed separately on 6 spectra whose likelihood are combined, as described in the text). Right panel: negative profile of the combined log likelihoods of 
RUN I and RUN II before (black) and after (color) the systematic uncertainty is included. 



In order evaluate the performances of the two approaches, 
we compared the distribution of the relevant statistical estima- 
tors for a thousand toy Monte Carlo simulated spectra for each 
CUORICINO detector, generated with rates corresponding to 
those directly measured in RUN II (see table 4) and no OvySyS sig- 
nal. Several pieces of information have been obtained from this 
comparison: 

• Both methods lead to compatible results. Figure 7 demon- 
strates a strong correlation between them. 

• Both are unbiased. Figure 8 (left panel) shows that for both 
methods, the distribution of the best fits divided by their 
statistical errors is compatible with a gaussian centered at 
zero with a variance equal to one. This is an important 
result which is not always guaranteed by maximum likeli- 
hood methods applied to low statistics. 

• Both have similar sensitivities. The distributions of the 
90% confidence intervals (figure 7 top panel) show a sen- 
sitivity of nearly 2.5x10^"* y, evaluated as the median of 
the distribution, following the Feldmann and Cousins pre- 
scription [28]. 

The last point offers a nice synthesis for a better understanding 
of the numbers we are dealing with. A wide range of limits can 
be reached in experiments with the same true background level 
(figure 8, right panel) and therefore with the same sensitivity. 
In this respect, quoting only the limit reached by an experiment 
can be misleading if the sensitivity is not also mentioned. 

Having such an impressive correspondence between the two 
approaches, the choice of the one or the other is somewhat ar- 
bitrary. We opted for the third method because it is consistent 
with the previous analysis of CUORICINO data [12, 30] and 
because of its intrinsic simplicity. 



7. Systematic uncertainties 

In our analysis we have identified the following sources of 
systematic uncertainties: 

• the calibration uncertainty; 

• uncertainty in the signal efficiency; 

• the background shape; 

• the energy window of the fit. 

In the case of the calibration uncertainty we have a direct esti- 
mate of its magnitude coming from a dedicated analysis of the 
residuals (see section 3): we reconstruct the position of a peak 
in the Qvfifi region with a precision AE - +0.8 keV. This sys- 
tematic error has been included directly in the fit as a gaussian 
fluctuation in the energy position. Since this uncertainty is sig- 
nificantly larger than the error (0.013 keV) on the Q-value [17], 
the systematic error discussed here also automatically includes 
the experimental uncertainty on the measured '^"Te transition 
energy. 

To evaluate the uncertainty from the choice of energy win- 
dow and shape of the background spectrum, we varied the 
model for the background (flat, linear and parabolic) at four dif- 
ferent energy intervals centered at the Q-value of the Ov/?/? de- 
cay (we increased the lower and the upper bound of an energy 
window starting from 2527±30 keV at increasing steps of 5 
keV). The results of this analysis have shown an average varia- 
tion inF"" of 3 X 10-' V' ■ 

The uncertainty on the signal efliciency is reported in table 2 
to be 1.1% for the big crystals and 1 .4% for the small crystals, 
both of which are negligible compared to the contributions from 
the energy scale and background parametrization uncertainties. 
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Table 4: Background and Co rates expressed in counts/keV/kg/y obtained by the combined fit of RUN I and RUN II data. The absolute rates are about 20% higher 
because the efficiencies have not been yet included. The high background rate obtained for '^"Te-eniiched detectors is mainly due to a higher intrinsic contamination 
of these crystals. 



Fit Parameter 



big crystals 



small crystals 



"TOT 



Te-enriched crystals 



RUN I: flat background rate 
RUN 11: flat background rate 
RUN 1: ^°Co rate 
RUN 11: '^°Co rate 



0.20±0.02 
0.153±0.006 
4.6+1.5 
2.5+0.3 



0.20±0.02 
0.17+0.02 

9±6 
1.7 + 0.8 



0.8±0.4 
0.35+0.05 
0± 14 
+ 3.5 



Table 5: We compare the most stiingent 90% C.L. half-life lower limits present in literature (column 2). For each experimental result (rows 1-4) and for each of the 
considered NME evaluations (column 3-6), we report a m„ range. This identifies the upper bound on the neutrino Majorana mass according to the different results 
reported by the same author (when varying some of the parameters in the used nuclear model). In the two lower rows we compare the nigg range obtained for the 
95% C.L. half-life limit on '''"Te (row 5) with the positive signal quoted by [14] (row 6). For this last case the range corresponds to the 2 sigma range in the 
measured half-life. 



Isotope 


' 1/2 

[y] 


QRPA [5] 
iriee [meV] 


QRPA [6] 
niee [meV] 


SM [7] 
niee [meV] 


IBM [8] 
[meV] 


i^OXe (CUORICINO, this work) 
^*Ge (Heidelberg-Moscow collaboration [9]) 
i™'Mo (NEMO coflaboration [11]) 
i^^Xe (Dama/LXe [13]) 


> 2.8 10^4 

> 1.910^5 

> 5.8 10^3 

> 1.2 1024 


< 300-570 

< 230^00 

< 610-1270 

< 700-1640 


< 360-580 

< 280-460 

< 810-1430 

< 800-1230 


< 570-710 

< 530-640 

< 1020-1270 


< 350-370 
<270 

< 830-850 

< 640-670 


i^"Te (CUORICINO, this work, 95% C.L.) 
^^Ge (Heidelberg-Moscow experiment [14]) 


> 2.3 lO^* 

- 9 9Q+O.881 n25 


< 340-630 
= 180^30 


< 390-640 
= 220-500 


< 620-780 
= 410-700 


< 390-410 
= 210-290 



8. The CUORICINO final result 

As a next step, we added the contributions from big, small 
and enriched crystals from RUN 1, combining their likelihoods 
with the RUN 11 data and using a similar reconstruction for the 
response function as described in section 6. 1 . 
The background rates are shown in table 4 while in figure 9 
(left panel) we show the best fit and the corresponding 68% and 
90% C.L. limits. Figure 9 (right panel) shows the logarithm 
of the combined likelihoods of RUN 1 and RUN 11 before and 
after the systematic uncertainties are included. 

The resulting best fit for the OvjSyS rate of '^"Te is: 

r{^^^^ = (-0.25 + 1.44(stat) ± 0.3(syst)) x lO^^Sy-i 

This result is compatible with zero, and the corresponding 
90% C.L. half Hfe lower bound is: 



r^;2 > 2.9 X 1024y 



90% C.L. limits that dififerent experiments (with the same expo- 
sure and sensitivity) can yield (figure 8). This is the reason we 
prefer to quote the sensitivity of the experiment together with 
the limit. It should to be mentioned that, in reference [12], we 
used an older value of the '^'^Te transition energy energy which 
had a much larger error and a slightly higher central value than 
the recently measured one [17, 31]. With the same data used 
in [12], using the new result for the transition energy (with its 
smaller error) pushes the limit toward a lower half-life. 

The inclusion of the systematic error modifies the likelihood 
profile for our data as shown in the right panel of figure 9. The 
profile can be considered the of our fit as a function of all 
the possible F"''. Thus we will refer to it as xltat- If we adopt 
the hypothesis that our knowledge of F"'' is smeared — near 
the best fit values — by a gaussian systematic uncertainty of 
magnitude (Tsyst, the total x^ot wiU be: 
1 1 1 



Xtol 



2 2 

X slat X sysl 



(8) 



where the simplest approximated form of x„st 



This limit is almost identical to the one we published in 
[12], despite the increase of the total exposure by a factor ~1.6. 
In general it is expected that the limit scales with the square 
root of the exposure, i.e. we would expect an improvement of 
about a factor 1.3 which is by far smaller than the spread in the 



X s\st 



(F--F«:j 



(9) 



With this modification of x^ and because the systematic un- 
certainty is small compared to the statistical error, we obtain a 
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slightly weaker limit on the half life: 

r'i'jj > 2.8 X lO^^y 

As it is a standard approach in OvySjS literature, we also present 
the 95% C.L. limit on T*j*j2 including systematic uncertainties: 

T^l% ^ 2.3 X lO^'^y 



9. Conclusion 

In this paper we have presented the CUORICINO final result 
on OvySyS in '-"'Te , obtained with an exposure of 19.75 kg-y of 
^■"'Te, including a detailed study of systematic errors for the first 
time. A half fife Hmit of 2.8 lO^^ y at 90% C.L. is obtained (2.9 
10^"^ y if systematic errors are not included), to be compared 
(as discussed in section 6.2) with an experimental sensitivity** 
of ~2.6 10^* y. This limit can be used to extract an upper limit 
on niec using the theoretical NME evaluation for '■'"Te nucleus. 
We report here results obtained using the most recent nuclear 
calculations found in literature: 

• 300-570 meV using the Quasiparticle Random Phase Ap- 
proximation (QRPA) evaluations of reference [5] 

• 360-580 meV using the QRPA evaluations of reference [6] 

• 570-710 meV using the Shell Model (SM) evaluations of 
reference [7] 

• 350-370 meV using the Interacting Boson Model (IBM) 
evaluations of reference [8] 

Note that, for each reference, a range (and not a single value) 
for nice is presented, reflecting the different results for the NME 
obtained by the authors when varying model parameters, such 
as the treatment of the short range correlations or the value of 
gA (the axial-vector coupling). Then, the interval 300-710 meV 
can be taken as the final range for the 90% C.L. upper bound on 
niee (at 95% C.L. this becomes 340-780 meV). 

In table 5 we compare this result with the most stringent 
90% C.L. half-life lower limits present in literature. For each 
experimental result we report the nig,, range obtained with the 
NME evaluations here considered. Despite the differences be- 
tween the NME evaluations, it is evident that CUORICINO is 
one of the most sensitive experiments performed to date. 

Finally the bottom two rows in table 5 compare the 95% 
C.L. half-life limit on '^"Te obtained in this work with the 2 
sigma range corresponding to the positive signal quoted by 
[14] and obtained with a re-analysis of the Heidelberg-Moscow 
data. The two results are clearly compatible. 



^This is the sensitivity evaluated for tfie total (RUN I + RUN II) statistics. 
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